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Abstract 



The cosmological QCD transition affects primordial density perturbations. If 
the QCD transition is first order, the sound speed vanishes during the transi- 
tion and density perturbations fall freely. For scales below the Hubble radius 
at the transition the primordial Harrison-Zel'dovich spectrum of density fluc- 
tuations develops large peaks and dips. These peaks grow with wave number 
for both the hadron-photon-lepton fluid and for cold dark matter. At the 
horizon scale the enhancement is small. This by itself does not lead to the 
formation of black holes at the QCD transition. The peaks in the hadron- 
photon-lepton fluid are wiped out during neutrino decoupling. For cold dark 
matter that is kinetically decoupled at the QCD transition (e.g., axions or 
primordial black holes) these peaks lead to the formation of CDM clumps of 
masses IO-^^Mq < Mdump < W'^^Mq. 
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I. INTRODUCTION AND RESULTS 



In the early universe a transition from a quark-gluon plasma (QGP) to a hadron gas 
(HG) took place at a certain transition temperature T^,. At high temperatures asymptotic 
freedom of QCD predicts the existence of a deconfined phase, the quark-gluon plasma. At 
low temperatures quarks and gluons are confined in hadrons. The deconfined phase at high 
temperatures was found by lattice QCD simulations The order of the QCD transition 
is still under debate, but recent lattice QCD results for three quark flavors with massive 
strange quarks indicate that the QCD transition is first order 0. The value of the latent 
heat is available only for quenched lattice QCD (gluons only, no quarks) These data 
give a small latent heat /; [l/T^) is about one fifth of the difference in entropies between an 
ideal massless quark-gluon plasma resp. hadron gas. Two-flavor calculations give ^ 140 
MeV 0]. Ongoing and future heavy ion experiments at the world's largest colliders (SPS, 
RHIC, LHC) are searching for the QGP §. 

At the QCD transition the Hubble radius i?H ~ n^-pi/T^ is about 10 km. The mass 
inside a Hubble volume is ~ IMq. The transition lasts about a tenth of a Hubble time, t^, 
with the latent heat from quenched lattice QCD. The Hubble time at the QCD transition 
is tn = Ru/c ~ 10~^ s. This is extremely long in comparison with the relaxation time scale 
of the strong interactions, which is about 10~^^ s. Thus, the transition is very close to an 
equilibrium process. 

For a flrst-order QCD transition hadronic bubbles nucleate during a short period of 
supercooling. The crucial parameters for supercooling are the surface tension a and the 
latent heat I. Data for a and / are available only for quenched QCD |^ and give a very 
small surface tension. This gives a very short period of supercooling, At^c ~ 10~^tH; and 
a typical bubble nucleation distance dnuci ~ 1 cm ^ 10~^-Rh for homogeneous nucleation 
I^J^. The hadronic bubbles grow very fast, within 10~Hii, until the released latent heat has 
reheated the universe to T^. For the remaining 99% of the transition the HG and the QGP 
coexist at the pressure Phg(^*) = Pqgp(^*)- During this time the hadronic bubbles grow 
slowly and the released latent heat keeps the temperature constant until the transition is 
completed. The energy density decreases continuously from pqgp(T^) at the beginning of 
the transition to pnciT^,) when the transition is completed. 

In the mid 80s interest in the cosmological QCD transition arose because it was realized 
that a strongly flrst-order QCD phase transition with large surface tension could lead to 
observable signatures today. Witten p pointed out that the separation of phases during 
the coexistence of the hadronic and the quark phase could gather most baryons in (strange) 
quark nuggets. However, the quark nuggets, while cooling, lose baryons, and evaporize, 
unless they contain much more than 10"^^ baryons initially 0. This number should be 
compared with the number of baryons inside a Hubble volume at the QCD transition, which 
is 10^°. Thus, the mean bubble nucleation distance should be > lO'^Ru ~ 100 m, which 
is about a factor 10"^ too big compared to the nucleation distance suggested using recent 
lattice results 

Applegate and Hogan found that a strongly flrst-order QCD phase transition with large 
surface tension could induce inhomogeneous nucleosynthesis fl^, assuming a nucleation 
distance which is ruled out by recent lattice results . The large isothermal baryon fluctua- 
tions, induced during the separation of phases, could lead to inhomogeneous initial conditions 
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for nucleosynthesis. As a minimal requirement for inhomogeneous nucleosynthesis the mean 
bubble nucleation distance should be larger than the proton diffusion length, which corre- 



sponds to ~ 3 m W^ at the QCD transition. This is more than 10 times larger than the 



nucleation distance based on lattice results . In addition, the observed cosmic abundances 
of light elements do not allow inhomogeneous nucleosynthesis, except in a small region in 



parameter space corresponding to an inhomogeneity scale of ~ 40 m [TI]. 

In a recent letter [|T^ and in the present work we look at matter averaged over scales A 
much larger than the bubble separation, e.g. A of order of the Hubble radius i?H- Therefore 
we deal with the bulk properties of the fluid in contrast to previous investigations which dealt 
with physics at the bubble separation scale. We found that the evolution of cosmological 
density perturbations is strongly affected by a first-order QCD transition for subhorizon 
scales, A < -Rh- Preexisting cosmological perturbations with an approximate Harrison- 
Zel'dovich spectrum are predicted by inflation |1T^-[17| and have been observed in the 
temperature fluctuations of the cosmic microwave background by the COBE satellite [0. 

We showed in that pressure gradients and the isentropic sound speed (for wavelengths 
A much larger than the bubble separation), = idp/dp)^'^ , must be zero during a first- 
order phase transition of a fluid with negligible chemical potential (i.e. no relevant conserved 
quantum number). The sound speed must be zero, because for such a fluid the pressure can 
only depend on the temperature, p{T), and because the transition temperature T^, has a 
given value, it cannot depend on any parameter, hence p(T-t) = is a given constant, and 
Cs = 0. This is in contrast e.g. to the water-vapor system, where the number of molecules is 
conserved, the pressure also depends on the particle density in a given phase, p{T,n), and 
where a two-phase system has an equilibrium pressure (the vapor pressure) which depends 
on the temperature, p*(T), hence Cs ^ 0. During the entire QCD transition the sound speed 
stays zero and suddenly rises back to the radiation value = 1/ a/3 after the transition 
is completed. Pressure varies continuously and goes below the ideal radiation fluid value 
p = p/3, but stays positive. Jedamzik ||19| independently pointed out that a significant 
decrease in the effective sound speed Cg during the cosmological QCD transition is expected. 

Since the sound speed is zero during a first-order QCD transition, there are no pressure 
gradients and no restoring forces. Preexisting cosmological perturbations go into free fall. 
The superhorizon modes (at the time of the transition), A > i?H, remain unaffected. The 
spectrum of subhorizon perturbations develops peaks in bpjp which grow with the wave 
number. The details of the enhancement depend on the QCD equation of state near T^,. 
The peaks grow at most linearly in wave number. The subhorizon peaks arise because 
preexisting acoustic oscillations at shorter wavelengths have a higher velocity (for a given 
5p/p) at the beginning of the transition and therefore get a larger amplification factor for 
during the transition. 

Photons and relativistic leptons, i.e. the three neutrinos, and /i^, also contribute to 
the energy density at the temperatures of interest. Weak and electro-magnetic interactions 
are fast enough to keep the photons and leptons in thermal equilibrium with the QGP/HG. 
Therefore the collective behavior of all these coupled particles can be described by one 
perfect (i.e. dissipationless) fluid, which we shall denote as radiation fluid (RAD). 

As a second fluid we include cold dark matter (CDM) which dominates the energy 
density of the universe today. At the time of the QCD transition p^°^(Tj ~ \^-^ p^^^iT^), 
therefore the gravity generated by CDM can be neglected. We must distinguish two types of 
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CDM with respect to the issue of kinetic decouphng at the time of the QCD transition. The 
hghtest supersymmetric particle, which is probably the neutralino is tightly coupled 
to the radiation fluid (and is treated as part of the RAD fluid) at T ~ T^. The scattering 
rate of the neutralino is of the order of the weak interaction scattering rate, the neutralino 
decouples kinetically at a much lower temperature of a few MeV. Kinetic decoupling must 
be distinguished from freeze out (when the pair production and annihilation rates fall below 
the Hubble rate), which happens way before the QCD transition. In contrast the second 
type of CDM is kinetically decoupled at the QCD transition, e.g. axions or preexisting 
primordial black holes. Note that in our figures and equations kinetically decoupled CDM 
will be labeled by CDM for brevity. 

Kinetically decoupled CDM falls into the potential wells provided by the dominant ra- 
diation fluid. Thus, the CDM spectrum is amplified on subhorizon scales. The peaks in 
the CDM spectrum go nonlinear shortly after equality. This leads to the formation of CDM 
clumps with mass < 10~^^ Mq. Especially the clumping of axions has important implications 
for the axion searches using the magnetic field induced axion decay ^1 . 



The formation of primordial black holes (PBHs) should be particularly efficient during 
the QCD epoch due to a substantial reduction of pressure forces during adiabatic collapse 
as pointed out by Jedamzik |TP[. The PBH mass function is expected to exhibit a pro- 



nounced peak on the QCD- horizon mass scale ~ IMq. He proposes that these black holes 



could account for massive compact halo objects (MACHOs) observed by microlensing P2 



However, we found by our linear analysis [|T^ that the amplification of fluctuations at the 
QCD horizon crossing scale is only a factor 1.5. For standard models of structure forma- 
tion without tilt, the amplitudes are not big enough to produce a cosmologically relevant 



amount of black holes ||23|. A tilted spectrum could be fine tuned to produce black holes at 
the QCD scale, but the spectrum would need a break just below the QCD scale in order not 
to overproduce smaller black holes. With the need of such a doubly fine-tuned preexisting 
spectrum the main effect would not be due to the QCD transition. 

The infiuence of the QCD transition on primordial gravitational waves has been inves- 
tigated by one of the present authors in Ref. It was found that the dramatic drop in 



relativistic degrees of freedom during the QCD transition (from 51.25 in the QGP to 17.25 
in the HG) induces a step by 30% in the primordial spectrum of gravitational waves. Today 
this step is located at ~ 10^^ Hz. Gravitational waves with that frequency could be detected 
by pulsar timing, however for gravitational waves from infiation the amplitudes are far too 
small to show up in the pulsar timing residuals. 

This paper is organized as follows: A discussion of hot QCD, bubble nucleation, and the 
vanishing of the sound speed for wavelengths A much larger than the bubble separation can 
be found in Sec. II. In this section we present a fit to recent lattice QCD results. We give 
the thermodynamical curves for w = p/p and the sound speed as a function of the scale 



factor at the QCD transition for three cases: The bag model, [^, which gives a simple 
parameterization, our fit to lattice QCD and a smooth crossover. 

In Sec. HI we derive the amplification of density perturbations at the QCD transition 
both numerically and analytically. We present the evolution equations for cosmological 
perturbations in a first-order form in uniform expansion gauge, which gives the same simple 
structure for the general relativistic dynamics as for the Jeans equations. The details of 
the amplification depend on the equation of state at the phase transition. We present three 
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cases: First, we use the bag model |2^, which allows a simple discussion of the effects. In 
the bag model, the amplification of subhorizon perturbations grows linearly in wave number 
k. Next, we use our fit to recent lattice QCD data p6| , p7| , which indicate a smaller latent 
heat P]. The amplifications for the lattice fit grow proportional to A;^/^. As a third example 
we study the possibility of a smooth crossover I^Hj. In this case, the sound speed stays finite 
and subhorizon scales are not amplified. The spectrum for a crossover has peaks only around 
the horizon scale. We prove that the spectra of primordial perturbations are not affected 
for scales greater than the Hubble horizon at the transition. 

In the last section (Sec. IV) we discuss observable consequences of the large peaks and 
dips in the spectrum of density fluctuations. Collisional damping around the time of neutrino 
decoupling wipes out all subhorizon inhomogeneities in the radiation fluid before big bang 
nucleosynthesis (BBN), and homogeneous BBN cannot be affected. We discuss black hole 
formation during the QCD transition and conclude that it is highly unlikely unless the 
primordial spectrum is fine tuned at the QCD scale. 

The most interesting prediction is the clumping of those types of CDM which are ki- 
netically decoupled at the QCD transition and do not suffer from collisional damping by 
neutrinos. E.g. axions, if they are the CDM, are not distributed smoothly within the halo 
of our galaxy, they come in clumps. 

Throughout the paper we assume that our universe is spatially fiat (i7 = 1) and most 
of its matter is CDM today. Whenever we give numbers we use a Hubble parameter of 
Ho = 50 km s"^Mpc~\ 

II. THE COSMOLOGICAL QCD TRANSITION 
A. Hot QCD 

The baryon number density in the universe is extremely small compared to the entropy 
density, i.e. hb/s = O{10~^^). At the time of the QCD transition the chemical potential for 
quarks resp. baryons is negligible as far as the pressure p and energy density p are concerned, 
p and p depend on the temperature only. 

The behavior of p{T) and p{T) near the QCD transition must be given by nonperturbative 
methods, lattice QCD. In Fig. |l] we plot lattice QCD data for p{T) and p{T) divided by p 
of the corresponding ideal gas. We show the lattice results for two systems, quenched QCD 
(no quarks) [^], and two- flavor QCD |2^. For quenched QCD the lattice continuum limit 
is shown. For two flavor QCD the data with six time steps {Nt = 6, a ~ 0.2 fm) and a 
quark mass arrig = 0.0125 is shown. This corresponds to a physical mass rUg ~ 14 MeV, a 
bit heavier than the physical masses of the up and down quarks. On the horizontal axis we 



plot (T/T^). We note that the critical temperature for quenched QCD is ~ 260 MeV [26 
and for two-flavor QCD Tj, ~ 140 MeV [Q. Unfortunately p{T) and p{T) for three quarks 
(u,d,s) with physical masses is not available yet. For T/T^ = 4 energy density and pressure 
for quenched QCD are still 10% resp. 15% below the ideal gas value. It is remarkable that 
p/pideai and p/pideai vcrsus T/T^ is quite similar for quenched QCD and two-flavor QCD. 
Moreover, the temperature dependence of the rescaled pressure for four-flavor QCD ||29| is 
quite similar to quenched QCD. 
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At temperatures below quarks and gluons are confined to hadrons, mostly pions. At 
present the hot pion phase is not seen in the two-flavor lattice QCD, since the pion comes 
out too heavy (0.3 < m^/m^ < 0.7 from ||^, whereas the physical ratio is 0.18). 

The second law of thermodynamics connects pressure and energy density, and reads for 
a fluid without chemical potential (no relevant conserved quantum numbers) 

P^r|-p. (2.1) 

The Maxwell relation for the free energy gives the entropy density, 

p..) 

From homogeneity the free energy density f{T) = —p{T), and f{T) contains the full ther- 
modynamic information. The lattice QCD results for s/sqq^p are shown in Fig. |^. 

The bag model gives a simple parameterization for p{T) and hence for p{T), s{T). The 
bag model ansatz represents the short- distance dynamics by an ideal gas of quarks and 
gluons and the long-distance confinement effects by a constant negative contribution to the 
pressure, the bag constant B, 

PqgAT) = p'S.iiT) - B . (2.3a) 

We include gluons and u,d quarks, which are effectively massless a.tT ^ T^, hence Pqgp(T) = 
(7r^/90)(7QQpT^, where g* is the effective number of relativistic helicity states, g* = S'bosons + 
7/8(?fermions, S'qgp = 37. The low temperature phase is a hadron gas (HG). We model it as 
an ideal gas of massless pions, 

PHG = ^^^HG^', (2-3b) 

because the contribution of pions is small anyway, 5'hg/5'qgp — 3/37. 

At the phase transition (T = T^) the quark-gluon phase and the hadron phase have the 
same pressure, p = p^,. The stable phase for T 7^ is the one which gives the minimum 
for the free energy density /(T), hence a maximum for p{T). The condition for pressure 
equilibrium at gives the relation between and B; using ( |2.3a] ) and ( [^.3b| ) we obtain 



B = ^{ghop - 9ug)T:. (2.4) 

We take from lattice QCD calculations with two flavors, which give G (140, 160) MeV 
|l|,^. This corresponds to a range of bag constants 5^/^ G (195, 221) MeV. This range is 
in agreement with fits to the light-hadron masses, which yield 5^^^ = 145 to 245 MeV (a 
compilation of various bag model light-hadron fits can be found in [0). We adopt the value 
T, = 150 MeV. 



Energy density and entropy density in the bag model follow from Eqs. ( p.l|) - (|2.3a|) 



Pqgp=PqgI + B (2.5a) 
^QGP = 4gp. (2-5b) 
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with p''^'"'^ = {n'^/30)g*T^ and s''^'"'^ = {271"^ / Ah)g*T^ for massless particles. The bag results 
for p/pqgp, p/Pqgp5 ■^/sqqp are shown in Fig. |^ and Fig. ||. Note that s/T^ is a simple step 
function in the bag model. The QCD transition is first order in the bag model, and the 
latent heat per unit volume / = T^As, 



IS 



45 



(2.6) 



where A^f = ^fqcp - guG- 

The order of the QCD phase transition and (for a first-order transition) the magnitude 
of the latent heat is still a subject of debate. In quenched QCD the phase transition is of 
first order The latent heat was determined to be Z ~ 1.4T^ It is useful to take the 
ratio Rl of the latent heat to the value T^.As^'^'^^^, where As^'^*^'^' is the difference in entropy 
between an ideal massless HG and an ideal massless QGP, 



Rl = 



I 



(T,As) 



ideal 



1 bag model 

0.2 quenched lattice QCD 



(2.7) 



For two light quarks it is likely that the transition is a crossover [| . This is in agreement 
with theoretical considerations which predict a second-order phase transition in the 
massless quark limit. For three light flavors and for the physical relevant case of two light 
and one heavy flavor the phase transition is likely to be of flrst order 0. This result was 
obtained using the Wilson quark action, whereas results with staggered quarks ||3^ indicate 
a crossover for the physical quark masses. For four quark flavors the transition is flrst order 
| 29| . Since the latent heat for lattice QCD is known for quenched QCD only, we decided 
to use the latent heat ratio Rl = 0.2 from quenched QCD as an indication for the physical 
case. 

We also need an analytic representation of the lattice QCD data. We decided to flt the 
entropy density (Fig. |^) and to derive the other thermodynamic quantities with Eqs. ( ^.11 ) 
and ( p.2| ). Below we consider again an ideal gas of the three pions. We flt the shape of 
the lattice QCD data by 



(27rV45)T3 



g*^^ + Ag*9{T-T^ 



Rr + a- Rl 



1 



T 



(2.8) 



A good flt for our purpose is obtained for 7 G (0.3, 0.4). We flx 7 = 1/3. The flt is plotted 
in Fig. |. 

Finally, we do not neglect the possibility that the QCD transition could be a smooth 
crossover. Such a crossover has been modeled by a simple interpolation between the ideal 
values of the entropy density with a tanh. It was flrst given in Ref. ||28| and is used in 
relativistic hydrodynamical simulations of heavy- ion collisions |^ . The entropy density is 
written as 



^crossover 



^ideal 
'QGP 



2 fi'QGP 



tanh 



AT 



(2.9) 



We choose AT/T^ = 0.1, see Fig. ^ This gives the shape of the initial rise of the lattice 
data, but it necessarily fails to give the slow rise at higher temperatures from s/sqq^p ~ 0.6 
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to 1. Comparing this model to the bag model and the lattice fit allows us to identify 
phenomena caused by the cosmological QCD transition that are independent of the order 
of the transition. The effects for wavelengths close to the Hubble scale will be shown to be 
approximately independent of the order of the transition. 

The early universe at T around 150 MeV also contains photons and relativistic leptons 
(e^,/i^, neutrinos) with a pressure 

Pi7 = ^9!,T' . (2.10) 

where = 14.25. We use the subscript minus to denote the value of some quantity at 
the beginning of the transition, i.e. when the universe has cooled to from above, and the 
subscript plus to denote the value of the same quantity at the end of the transition, i.e. 
when the temperature starts to decrease again. 



B. Bubble nucleation 

The expansion of the Universe is very slow compared to the strong, electro-magnetic, 
and weak interactions around T^. To be more explicit, the rate of the weak interactions is 
Tw ~ G'|Tf ^ 10~^^ GeV, the rate of the electro-magnetic interactions is Fem ~ a^^^r ~ 10^^ 
GeV, and the rate of the strong interactions is Fs ~ (T^)T^ ^ 10"^ GeV. These rates have 
to be compared to the Hubble rate H ~ Tf/mpi ~ 10~^^ GeV. Thus, leptons, photons, 
and the QGP/HG are in thermal and chemical equilibrium at cosmological time scales. All 
components have the same temperature locally, i.e. smeared over scales A ~ 10~^i?H- At 
scales A > 10^'' -Rh strongly, weakly, and electro-magnetically interacting matter makes up 
a single perfect (i.e. dissipationless) radiation fluid. 

In a first-order phase transition the QGP supercools until hadronic bubbles are formed 
at < P, |5^J53| . Without 'dirt' the bubbles nucleate due to thermal fluctuations 
(homogeneous nucleation). The probability to nucleate a bubble by a thermal fluctuation 
is proportional to exp(AS'), where AS" is the change in entropy by creating a bubble. The 
second law relates AS* to the minimal work done in this process, which is the change in the 
free energy because the volume and temperature are fixed [^]. The change in free energy 
of the system by creating a spherical bubble with radius R is 

4:71 

AF = —{pQGP - Phg)R^ + 47ra/22 , (2.11) 

where a is the surface tension. Bubbles can grow if they are created with radii greater 
than the critical bubble radius -Rcrit- Smaller bubbles disappear again, because the free 
energy gained from the bulk of the bubble is more than compensated by the surface energy 
in the bubble wall. -Rcrit is determined from the maximum value of AF{R), Rcrit{T) = 
2cr[pHG(^) ~ PQGpiT)]~^ . At the critical bubble size diverges, and no bubble can be 
formed. Finally, the probability to form a hadronic bubble with critical radius per unit 
volume and unit time is given by 

/(T) = /oexp(-^) , (2.12) 
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with AFcrit = 167r(T^/[3(j9HG — Pqgp)^]- For small supercooling i] = 1 — T/T^, ^ 1 we may 
evaluate {puG ~Pqgp)(T) by using the second law of thermodynamics, i.e. phg ~Pqgp ~ 
and thus 

AFe« « — ^ . (2.13) 

For dimensional reasons the prefactor Jq ~ T^. A more detailed calculation of Jq within the 
bag model has been provided in It was shown in Ref. [0 that the temperature depen- 
dence of the prefactor Iq can be neglected for the calculation of the supercooling temperature 
Tsc in the cosmological QCD transition. Furthermore a purely numerical prefactor to 
would be irrelevant for our purposes as explained after Eq. ( p.l5| ) . Therefore the probability 
to form a critical bubble per unit volume and unit time can be written as 

I^T!exp(-A/r]'), (2.14) 



with A = lQ7va^/{3PT^). 

The surface tension o" is a crucial parameter for Tgc and -Rcrit- The absence of surface 
excitations in hadronic spectra suggests that a^^^ <^ 5^/^ [^. In lattice QCD rather small 
values are found. Refs. find that a ~ O.OlSTf for quenched lattice QCD. There are no 
values for unquenched QCD available yet. However, an upper bound was obtained for the 
case of four-flavor lattice QCD in Ref. [^, i.e., cr < O.IT^. Using the results from quenched 
lattice QCD we find A = 2.9 x 10"^ 

After the first bubbles have been nucleated, they grow most probably by weak de- 
flagration p^ , |39| , ^ ,|6|] . The deflagration front (the bubble wall) moves with the velocity 



'^defl l/V^ [0. The energy that is released from the bubbles is distributed into the 
surrounding QGP by a supersonic shock wave and by neutrino radiation. This reheats the 
QGP to and prohibits further bubble formation. The shock front propagates with a speed 
slightly above the sound speed. The amplitude of the shock is very small p3. Neutrinos 



have a mean free path of 10~''Th at T^. If they do most of the heat transport, it goes with 
^^heat = C'(c). A detailed calculation of the supercooling temperature in the cosmological 
QCD transition was given in Refs. [p5|,|6|. These calculations assume small supercooling. 



The suppression of bubble nucleation due to already existing bubbles is neglected. 

The supercooling temperature fraction rjsc turns out to be about the same for the 
schematic case of one single bubble nucleated per Hubble volume per Hubble time. 



A 



1/2 

^4x10"^ (2.15) 



or the realistic case obtained below of one bubble nucleated per cm^ per 10~^ of a Hubble 
time, which needs a supercooling 20% larger. This small increase of 20% in r/sc makes the 
bubble nucleation rate larger by a factor 10^^. 

The time needed for the supercooling is given by Atsc/^H = ''7sc/(3c^) = C(10~^). The 
critical size of the bubbles created at the supercoohng temperature is 

i?crit(r/sc) ~ ^ ~ 30 fm. (2.16) 
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Bubbles present at a given time have been nucleated typically during the preceding time 
interval Atnuci = -^/(d//dt) = O(10~^). Using the relation between time and supercooling 
T], drj/dt = 3Cg/tH, we find At^^ci/tu = V^/i^^^^D- During this time interval each bubble 
has distributed released latent heat over a distance ^ Atnuci'^heat- This distance has a weak 
dependence on the precise value of r/sc, but the bubble nucleation rate increases exceedingly 
strongly with rj until one bubble per volume ~ (Atnuci'^hcat)^ is nucleated. Therefore the 
mean bubble separation is 

3 

rfnucl ~ t^heatAtnucl ~ ^^^H = ^(lO'^i^H) = O(lcm), (2.17) 

6Cg A 

where we used fheat = 0{c),Cs = 0{c), which gives a typical value for the nucleation dis- 
tance. In [0 the mean bubble nucleation distance due to thermal fiuctuations (homogeneous 
nucleation) is calculated to be less than 2 cm, whereas for inhomogeneous nucleation (the 
first bubbles form at impurities like topological defects or primordial black holes) it might 
be a few m. 

The above estimate of bubble separation applies if the limiting factor for quenching is the 
distribution of released latent heat by means of sound waves and by neutrino free streaming. 
On the other hand the limiting factor could be given by the rate of release of latent heat, 
i.e. by the bubble wall velocity v^cfi- Since the period of supercooling lasts about 1% of 
the time needed for completing the entire first-order phase transition, 1% of the QGP must 
be converted to HG in the process of sudden reheating to T^, hence the bubble radius at 
quenching must reach 0.2 of the bubble separation, -Rbubbie ~ 0.2(inuci- With -Rbubbie ~ 
t;dofiAtnuci and using the above relation d^^^^i ~ ^^heatAtnuci, we require Vdcn > 0.2t;heat for 
consistency. If v^eR is smaller than this, the limiting factor for quenching is the rate of 
release of latent heat by bubble growth, and the bubble separation is 

3 

<ucl ^ t^deflAtnucl ~ ^^Rh, (2.18) 

6Cg A 

i.e. the bubble separation will be smaller than the estimate Eq. ( p.l7| ). 

To summarize, the scales on which non-equilibrium phenomena occur are given by the 
mean bubble separation, which is about 10~^-Rh- The entropy production is tiny, i.e. 
AS/S ~ 10~^, since the supercooling is small ~ 10~^. After supercooling, which lasts 
lO'^tH, the universe reheats in Atnuci ~ 10 ^tn- After reheating the thermodynamic vari- 
ables follow their equilibrium values and bubbles grow due to the expansion of the Universe 
only. This is in striking contrast to bubble nucleation at the GUT transition, where the 
particle interaction rate is of the order of the Hubble rate and the bubble nucleation rate 
per unit volume is very much smaller than H^. 

C. Adiabatic phase conversion 

We now discuss the evolution of the temperature and other thermodynamic variables as 
a function of the scale factor a. The evolution of w = p/p for the bag model, our lattice 
QCD fit, and for the smooth crossover ( |2.9|) is shown in Fig. |^. 
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Entropy is conserved, apart from the very short stage of reheating (~ 10~^tH) after the 
first bubbles have been nucleated. This allows us to calculate T = T{a) from d [s(T)a^] = 0, 
i.e. 

^ = (2 19) 

din a ds/dT ' ^ ^ 

except for T = in the case of a first-order phase transition. In the bag model s ~ and 
therefore T oc 1/a for T ^ T^. In the case of the lattice fit and the crossover the entropy 
is given by Eq. (|2.8| ) resp. Eq. ( p.9|) plus the contribution of photons and leptons. Pressure 
as a function of the scale factor is obtained from s = dp/dT, using Eq. ( |2.19| ). The energy 
density, hence w, follows from Eq. (p.l|) . 

While the QGP and HG coexist in a first-order QCD transition the expansion factor is 
determined by entropy conservation. 



a_ 




(2.20) 



where the index — (+) denotes the value of a quantity at the beginning (end) of the coexis- 
tence epoch. In the bag model the Universe expands by a factor a+/a_ ~ 1.4 until all QGP 
has been converted into the HG, whereas for our lattice QCD fit ( p.8|) the Universe expands 
by a factor a+/a„ ^ 1.1. The growth of the scale factor is related to a lapse in cosmic time 
by din a = Hdt. In terms of Hubble time the transition lasts O.Stu for the bag model and 
O.ltH for our lattice QCD fit. 

During a first-order QCD transition, i.e. T = T^, the pressure p(T^) = is constant. For 
any first-order QCD phase transition p{a) is obtained from the first law of thermodynamics 
dp = — 3(p + p^,)da/a. The result for a_ < a < a+ is: 



p(a) = (p(a„)+p.)( — ) -P*, (2.21) 



wia) 




1 



-1 



(2.22) 



where w = p/p. After the end of the phase transition we have made the approximation 
of noninteracting and massless pions, leptons and photons, p = p/3, and therefore = 
p(a+)/3. 



D. The sound speed 

As explained in the introduction, the isentropic sound speed (for wavelengths A much 
larger than the bubble separation), given by 

must be zero during a first-order phase transition for a fluid with negligible chemical potential 
(i.e. no relevant conserved quantum number). This behavior is shown in Fig. ^ for the lattice 
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QCD fit and the bag model. In the bag model = 1/3 before and after the transition. For 
a crossover the sound speed does not drop to zero. 

A strong decrease in the sound speed is observed in lattice QCD for Nj = |^ and for 
Nf = 2 1^. From our lattice QCD fit we can calculate the sound speed for T > T^. From 
Eq. ( |2.1| ) and Eq. ( p.2|) the sound speed is given by = (dlns/dlnT)~^. Inserting (p. 81), 
we obtain 

oc (l - . (2.24) 

For an analytic discussion of the behavior of acoustic modes during the phase transition it 
is important to know as a function of conformal time rj = J a~^{t)dt. In the radiation 
dominated regime a oc t]. Up to small corrections, this holds true during the transition. 
Let us denote the value of conformal time at which the sound speed vanishes by ri_. From 
Eq. ( |2.19 ) we find (//_ — ri)/rj- oc (1 — T^^/Ty, which implies 



ci = C'[iv--v)/V-]~ , (2.25) 

with C2 = {3^-^[^; + A^*i?L]/[A^*(l - i?L)]}i/^/7. For 7 = 1/3 the sound speed goes to 
zero linearly in conformal time. 

We now discuss why the isentropic condition applies during the part of the phase tran- 
sition after the initial supercooling, bubble nucleation, and sudden reheating to T^. During 
the second part of the transition, which takes about 99% of the transition time, the fluid is 
extremely close to thermal equilibrium, because the time to reach equilibrium is very much 
shorter than a Hubble time, i.e. the fluid makes a reversible transformation. This can be 
seen as follows. Across the bubble walls local pressure equilibrium is established immedi- 
ately, J9QGP = Prg locally. Local temperature equilibrium, Tqgp = Thg is established by 
neutrinos, which have a mean free path of 10~''i?H, enormously larger than the bubble wall 
thickness, and a collision time much shorter than the Hubble time. This local pressure and 
temperature equilibrium can only be satisfied if p = and T = at the bubble walls. Over 
distance scales of the order of the bubble separation (~ 1 cm) pressure (and therefore also 
temperature) is equalized with the velocity of sound, and thereby the released latent heat 
is distributed. This pressure equalization is very fast compared to the Hubble expansion 
velocity wh ~ 10~®c at the 1 cm scale. When analyzing cosmological perturbations we shall 
consider wavelengths A > IO^'^-Rh, for which neutrinos are tightly coupled, T,y/k ^ 1, see 
Sec. IV. For these wavelengths the radiation fluid behaves as a perfect (i.e. dissipationless) 
fluid, entropy in a comoving volume is conserved, one has a reversible process. On the other 
hand, below the neutrino diffusion scale, A < 10~^Rh, acoustic oscillations are damped away 
before the QCD transition. 

We have been criticized in ^2| for the use of the isentropic condition. In |^ it was 



claimed that the isentropic condition leads to an infinite sound speed in the bag model. 
This claim was wrong, because it was based on non-relativistic hydrodynamical equations 
which neglect the pressure contribution to the momentum density and neglect that sound 



waves in a radiation fluid move with relativistic speed. In a revised version of their work ||43 
they correct these mistakes. However they still maintain that the isentropic condition does 
not hold, but that instead the volume fraction of the two phases is frozen. They do not give 
any argument why such a freezing should occur, and their freezing assumption contradicts 
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the rapid equilibration of pressure and temperature between the two phases explained in 
the previous paragraph. Their freezing assumption would also entail supercooling from the 
Hubble expansion in the presence of a finite volume fraction of bubbles. This is in dramatic 
conflict with the extremely rapid reheating to (within 10~^tH) once enough bubbles (~ 1% 
of volume) have been formed as discussed in Sec. IIB. 



III. PEAKS AND DIPS IN THE DENSITY SPECTRUM 

A. Evolution equations for cosmological perturbations in uniform expansion gauge 

The evolution of cosmological perturbations at the QCD transition is analyzed for perfect 
fluids. We linearize the Einstein equations, e.g. p(x, t) =: po{t) +6p{'x,t). It is an excellent 
approximation at early times to take the Friedmann- Robertson- Walker metric fiat, (is^ = 
—dt^ + a{tyd:K ■ rfx. The evolution of the background quantities is given by the Friedmann 
equations, 

H = -^Po (3.1) 
H= _4vrG(po+po), 

where the subscript denotes background quantities. 

The perturbed metric is decomposed in a scalar, a vector, and a tensor part. We are 
only interested in density perturbations, which means the scalar (longitudinal) sector. For 
a time-orthogonal foliation of space-time the perturbed metric can be written as 

ds^ = - (1 + 2a) dt^ + a{tf [5ij{l + 2(^) + 2(9^(9^7] dx'dx^, (3.2) 

with di = We follow Bardeen 1989 and focus on the geometric properties of the 
constant-t hypersurfaces (slices) S. The perturbation of the lapse function between the S's 
is given by a(x, t). The intrinsic curvature of E is generated by <y9(x, t). 



(3) 



R[T] = -4 A yp. (3.3) 



The extrinsic curvature ^ of S has 3 terms, the unperturbed term, the perturbation of the 
isotropic part, 5[tri^'*j(E)] = k, and the anisotropic part (= shear of normals on S), which 
is generated by x(x, t), 

K^s] = + - {d^d, - ^ A 5^)x, (3.4) 

where 9* = a~'^5^^di. These geometric properties of S are computed from the metric pertur- 
bations as follows, X = '^^7; 1^ = ""3 {ip — Ha) — Ax- Vanishing shear of ?i(S) means x = 
and vice versa, but note that x (resp. 7) contributes both to the anisotropic part and to 
the isotropic part of K^j. 

We choose a slicing of space-time with unperturbed mean extrinsic curvature, 
5[tri^'*j(S)] = K = 0. This implies that our fundamental observers, which are defined to 
be at rest on the slice S, tt(obs) = n{T,), have relative velocities, which in the mean over 
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all directions follow an unperturbed Hubble flow. If the coordinate choice (gauge choice) 
is such that the time coordinate t is constant on the slices S, the gauge is flxed to be the 
Uniform Expansion Gauge (UEG) = Uniform Hubble Gauge (UHG) |I6| . 

In the literature the most common gauges are the synchronous gauge, where the funda- 
mental observers are moving inertially (i.e. a = 0), the longitudinal gauge (= zero shear 
gauge, X = 0), where the fundamental observers expand isotropically, and the comoving 
gauge, where the fundamental observers sit in the fluid rest frame, M(obs) = M(matter). The 
uniform expansion gauge (UEG) is free of spurious gauge modes (in contrast to the syn- 
chronous gauge), it is nonsingular in the superhorizon limit and therefore easy to handle in 
numerical calculations (in contrast to zero-shear variables, as pointed out by Bardeen [0), 
and it has the two basic fluid variables Sp and f peculiar (in contrast to the comoving gauge). 

As fundamental evolution equations in UEG for perfect fluids we take the energy- 
momentum conservation V^T'^ = 0, i.e., the continuity equation and (in the longitudinal 
sector) the 3-divergence of the Euler equation of general relativity. Each fluid separately 
must fulflU 

dte = -3H{e + 'K)-Atlj-3H{po+po)a (3.5) 
dt4> = -SH^p - vr - (po + Po)(y , (3.6) 

where e = 6p = ST'^q, n = 6p = l/3(5[trT*j], Viip = 5T^^ for each fluid separately. Note that 
Viip = Si = momentum density (Poynting vector). The system of dynamical equations is 
closed by Einstein's i?°o"6quation, 

(A + 3if)a = 47rG(e + 37r) , (3.7) 

where the sum over all fluids is taken on the right-hand side. In addition we need (p/p) 
and (vr/e)s = c^. Figs. ^ and |^, for adiabatic perturbations. Equations ( |3.5D - ( |3.7| ) deflne 
our general relativistic Cauchy problem (linear perturbations, longitudinal sector, perfect 
fluids). The initial data (e, ■?/') can be chosen freely on Sue- 

Einstein's i?°Q-equation in linear perturbation theory in UEG is an elliptic equation. 
Its Green function is a Yukawa potential with exponential cutoff at the Hubble radius, 
a ~ exp(— pr), where /i^ = —3H, instead of Newton's 1/r-potential in the Jeans equa- 
tions. R^Q gives the relative gravitational acceleration (geodesic deviation), averaged over all 
directions, of two test particles, which are moving nonrelativistically in the coordinate sys- 
tem at hand. In UEG SR^q is determined by a, the perturbation of the lapse. In Newtonian 
physics the relative acceleration averaged over directions is given by A$Newton, i-e. R^q and 
A<l>Newton auswcr the same physics question. In UEG a plays the role of the gravitational 
potential. 

Our basic equations have exactly the same structure as the Jeans equations: two flrst- 
order time-evolution equations for each fluid, Eqs. (|3.5| ) and ( p.6|) , supplemented by one 
elliptic equation, Eq. ( p.7| ). The UEG is the only gauge with this structure and with Ein- 
stein's equation appearing only in the elliptic equation. In the subhorizon limit (drop the 
last term in ( |3.5| ) and the H-teim in (|3.7|)) and with p <^ p these equations automatically 
reduce to the Jeans equations. The UEG is singled out for two reasons: First, the relative 
velocity of the fundamental observers in the mean over all directions is an unperturbed 
Hubble flow both in the UEG of general relativity and in the Jeans analysis. Second, only 
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the mean over all directions appears in the continuity equation, in the 3-divergence of the 
Euler equation, and in the i?°Q-equation. 

Out of the geometric quantities (k, x, V^, «) only the lapse a appears in the equations 
needed to solve the dynamics in UEG. The spatial curvature ((^) and the shear of the 
fundamental observers (x) are automatically absent from the above equations for perfect 
fluids, since these geometrical properties of S are irrelevant for the dynamical question at 
hand, (v?, x) are not needed at all to solve the dynamics. One can compute and x at any 
time from the energy and the momentum constraints, i.e. from the G^q equation and the 
longitudinal contribution to the equation. 



A v? 
Ax 



(3.8) 
(3.9) 



This sequence, solving the dynamics before using the G^q and G'^^ constraint equations, is 
different from usual solution strategies in linear perturbation theory. 

Since the evolution equations are linear, and the background is spatially homogeneous, 
each spatial Fourier mode with comoving wavelength k := kP^^^a evolves independently. 
It is convenient to rewrite equations ( p.5| ) - ( p.7| ) in the dimensionless variables b ^ tj p 
(density contrast), = kp\ijsilj/ p. The variable ip is related to the fluid velocity (peculiar 
velocity) fpcc = ^^^l"^! and has the same order of magnitude as 5 on subhorizon scales, 
^phys ^ _ ^-1 jj^ UEG the system of evolution equations, written in terms of conformal 
time 0' = dr, = adt, now reads 



w)S 



3(1 + w)a 



n 



V^' + (1 - 3w)^ 



—c 



n 



l + w) 



n 



-a 



+ 1(1 + w) 



a 



-^(1 + 3^5, 



(3.10) 
(3.11) 

(3.12) 



with H = {da/dri)/a and w = p/p. The Friedmann equation and the continuity equa- 
tion for the background were used. The continuity and Euler equation refer to each fluid 
separately. In the 
cH(E/vr(/))/(E 



i^Vequation w = (E/P^^V(E/ P^^^), S = (E/ e(^^)/(E/ P^^^) and 



/ 



:(/) 



In our numerical analysis we have used the exact general relativistic equations ( |3.10 ) - 
( p.l2| ), but it is instructive to look at the superhorizon A^^^^^ ^ H^^ and at the subhorizon 
\p^ys ^ limit of equations ( |3.10| ) - ( p.l2| ) for general pressure and sound speed. 

For superhorizon perturbations with arbitrary equation of state (e.g. QCD transition), 
we only have to keep the lowest order in k/Ti. in the evolution equations. The Equations 
( p.lOj ) - (|3.12|) for a one fluid model simplify to 



1.' 

n 



[l + 3w)S 



n 



ip' + {1 - 3w)ip 



k J 

1 k 

3H 



5, 



where the superhorizon limit of the Poisson equation (1 + w)a = — 1( 
Note that the sound speed drops out in the superhorizon equations. 



(3.13) 
(3.14) 

l + 3c^)5 was inserted. 
For the growing mode 
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the terms on the left hand side of the Euler equation add and if) 1 5 (x k/H. Therefore, 

L J grow 

the right hand side of the continuity equation is neghgible and the solution for the growing 
mode of the density contrast is 

6^{k/7if = {k^^yjHf (3.15) 

for superhorizon perturbations in the uniform expansion gauge. This result is independent of 
the thermodynamic behavior. It says that the evolution of the growing superhorizon mode, 
if expressed in terms oi Ti (as opposed to expressing it in terms of t or a), is unaffected 
by a phase transition. Therefore the QCD-transition cannot affect superhorizon scales, and 
the spectrum of density perturbations stays flat for A'''^^^ ^ H~^. This result can also be 
formulated differently: using the energy-constraint ( |3.^ ) we find ip = const for the dominant 
mode on superhorizon scales in UEG. This is the 'conservation law' [|l^ for the quantity 
( = if + ^i^^^^ S for w —1, since ( ^ ip in the superhorizon limit. 

For subhorizon analysis, A^*^^^ ^ H^^, with arbitrary equation of state, we can drop H in 
the general relativistic Poisson equation ( p.7|) since \H\ = 0{H^) <^ ^phys- the continuity 




equation ( |3.5| ) the time dilation term (last term) can be omitted since it follows from the 
Poisson equation that a = 0{H'^/k'^^yJ6 <C 6. In the subhorizon limit Eqs. ( |3.10|) - (|3.12|) 
simplify to 

1 u 

(3.16) 

-il + w)^a (3.17) 
l{l + 3cl)6. (3.18) 

In the case of the QCD transition, the transition time is short compared to the Hubble 
time (t+ — t-) < tn = and therefore self-gravity of the dominant radiation fluid can be 
neglected in the Euler equation ( p.l7|) for this component. On the other hand, CDM moves 
in the external gravitational potential provided by the radiation fluid, and is coupled to the 
radiation fluid via the gravity term in the Euler equation. Therefore we have to keep the 
gravity term in the CDM Euler equation. 

For a purely radiation dominated universe w = = |, the solution for the density 
contrast S and the peculiar velocity i/j can be written in terms of spherical Bessel and 
Neumann functions [Q, 

S = xf{x) + ^—x'^f{x) (3.19) 
a;^ + 6 ax 

^ = /(^)' (3-20) 

-|- 6 ax 

where x = kp^ys/H = k/H. For a radiation dominated universe x = kr], f{x) = 
Aji{x/ ^/S) + Bni(x / /VS. The normalization is chosen in such a way that the sub- 
horizon modes, A^*^^^ ^ H^^, of 5 and VSip are acoustic oscillations with constant amplitude 
y/WTW. We will fix our initial conditions to be the growing mode in the radiation domi- 
nated superhorizon regime. 



16 



6 = A-x^ (3.21) 
6 

^ = A^x\ (3.22) 
54 



B. Numerical results 

The transfer functions for the radiation fluid and the cold dark matter (CDM) are calcu- 
lated by integrating the exact general relativistic set of equations (|3.10| )-( p.l2| ) throughout 



the QCD transition. The initial conditions are fixed at T = 100 GeV, where all scales under 
considerations were far above the horizon. We fix the initial amplitude of the radiation fluid, 
Ain, by the growing mode of an exactly radiation dominated solution, Eqs. ( p.21| ) and ( p.22D . 



The normalization constant is given by the COBE normalized Harrison-Zel'dovich spec 



trum |T^, but it drops out in the transfer functions. The initial conditions for CDM are 



obtained assuming adiabatic perturbations, i.e. the entropy per cold particle is unperturbed 

,RAD/^CDM) ^ 0. Since pCDM ^ ^^CDM ^ 1/^^^.^^ and PRAD OC (1/Komoving)'/', the 

adiabatic initial conditions for (J^dm written 



^ODM ^ S^RAD (3 23) 

The initial fluid velocities are equal. Since we are in the linear regime, each spatial Fourier 
mode evolves independently. 

For a first-order phase transition it is convenient to use the volume of a comoving box 
OC instead of time as the evolution variable. Thus, the derivatives in Eqs. (|3.10|) and 
( p.ll| ) have to be replaced, 1-L~^ {d / dif) = a{d/da). In sections II. C. and II. D. we have given 
w and in terms of the scale factor. In the case of a crossover we use the temperature T as 
the independent evolution variable, since temperature is decreasing strictly monotonically in 
time. The derivatives in Eqs. ( p.lO|) and (p.ll|) are replaced by H~^{d/dri) = —3cj.T{d/dT). 



We obtain p{T) by integrating the entropy density s(T) , Eq. 

The transfer functions are evaluated by evolving each mode from the initial temperature 
Tin = 100 GeV to the final temperature Tf = T^/10 = 15 MeV, i.e. until the peculiar 
velocity of CDM has redshifted away. At Tf the universe is exactly radiation dominated. 
The amplitude of the acoustic oscillations is A^^^ = (5rad + 3^rad)^''^- The final amplitude 
of CDM is ^CDM _ |^CDM|_ rjnj^g transfer functions for the radiation and the CDM fluid are 
defined by 

r^AD(/t, in ^ f) ^ A^^^{k)\T, ^3 24) 

Am 

T™(A;, in ^ f) ^ ^ /^^'^^ , (3.25) 

^in 

where A^ = A^^^. 

For the bag model the transfer functions are shown in Fig. ^. Both transfer functions 
show huge peaks on small scales. The different scales k = 27i/X are represented by the 
invariant CDM mass contained in a sphere with radius A/2, 
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^CDM^;^phys) ^ ^^CDM ^^^^ ^ ^3_2g) 

assuming that ^cdm ~ 1 today. The largest scales in Fig. ^ correspond to the horizon at 
Tf = T^/10. The CDM curve also shows the logarithmic growth of subhorizon scales of CDM 
in a radiation dominated universe. The CDM peaks lie on top of this logarithmic curve. 

The peak structure starts at a scale ~ 10~^Mq in CDM mass. This scale corresponds to 
the horizon scale at the QCD transition, Mh = ^p'"°^(-Rh/2)'^. The radiation energy inside 
the horizon at is ~ IMq, but it is redshifted as M^^^(a) ~ (aequaiity/a)M'^°^. Scales 
which are above the horizon at the QCD transition are not affected in accordance with the 
general proof in the previous chapter, Eq. ( p.l5|) . For scales below M^^^ ^ 9 x 10~^Mq the 



radiation peaks grow linearly in wave number. This linear growth comes from the fact that 
the vanishing sound speed during the QCD transition implies a vanishing restoring force 
in the acoustic oscillations on subhorizon scales. Therefore, the radiation fluid falls freely 
during the transition, with a constant velocity given at the beginning of the transition. The 
density contrast S^^^ grows linearly in time with a slope k. CDM is moving in an external 
potential provided by the dominant radiation fluid, and is pushed by the strong increase in 
the gravitational force during the transition. The highest peaks have k/ki ~ lO'^, because 
on smaller scales the acoustic oscillations are damped away by neutrino diffusion already 
before the QCD transition (see Sec. IV). 

For our fit to lattice QCD the transfer functions 7"^^°, T^^^ are shown in Fig. |^. The 
peak structure starts at the horizon scale Mi, but the asymptotic growth on small scales is 
different. It grows oc k^^^ starting at a scale ^2 which corresponds to an invariant CDM mass 
of M2 = 2 X 10~^^Mq. The asymptotic envelope for small scales is indicated by the straight 
line in Fig. |^. The highest peaks show an amplification of a factor 400 and correspond to 
wave numbers k/k2 ^ 3 x 10^. The smaller amplification compared to the bag model is due 
to the decrease of the radiation fluid velocity before the transition. This happens because 
the sound speed drops below 1/a/3 already before the transition. The values for M2 and the 
exponent of asymptotic growth depend on our specific choice of the latent heat {Rl = 0.2) 
and the thermodynamic behavior at T > T^, parameterized by 7 (7 = 1/3). In the next 
chapter we will show, that the exponent of the asymptotic envelope for 7 G (0, 1) varies 
between 1/2 and 1. 

The processed spectrum for a crossover. Fig. |^, shows a similar behavior as for the bag 
model and the lattice fit on superhorizon and horizon scales. The peak structure starts 
at Ml, but on subhorizon scales there are no peaks. The level of the subhorizon transfer 
function for the radiation fluid is reduced to 0.83. This comes from the damping of the 
acoustic oscillations during the time with 7^ 1/3. 



C. Analytic solution for the radiation fluid 

The time evolution for subhorizon modes, X^^^^ -c -Rh = H^^, at the transition can be 
solved analytically. For the dynamics of the radiation fluid (QCD, photons, leptons) one 
can neglect cold dark matter (CDM), since pCDM^^RAD _ ~ 10^^. The transition 

time is short compared to the Hubble time at the transition, (t+ — t_) < tn = H~^. For 
subhorizon modes we can neglect gravity during the whole transition as we will show below. 
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The damping terms in the continuity equation ( p.l6| ) and Euler equation (|3.17|) are absent 



in the purely radiation dominated regime. During the transition the damping terms can be 
neglected in view of the huge amplification for a first-order phase transition. Eqs. ( p.l6| ) 



and (|3.17| ) now read: 



§> -kip = (3.27) 
^' + clk6 = . 

Written as a second order differential equation for 6, this is just an oscillator equation 

S" + uj^6 = , (3.28) 
with the time-dependent frequency uj{ri) = kcs{ri). 

1. Bag model 

We first discuss the origin of large amplifications for small scales in the bag model. Before 
and after the QCD transition the radiation fiuid makes standing acoustic oscillations in each 
mode k with sound speed = 1/3 and with amplitude and Aout for the incoming and 
outgoing solution respectively, see Fig. H. The incoming solution for the density contrast 8 
and the peculiar velocity v^V' reads 

5 = -AinCos(u;(r7 - ?7_) + v9_) (3.29) 
= AinSin {uj{ri - r]^) + (y9_) . 

This solution is valid until the beginning of the transition at r/ = and denotes the 
phase of the oscillation at During the transition the sound speed is zero. There are 
no restoring forces from pressure gradients and the radiation fiuid falls freely. Since the 
duration of the transition is short compared to the Hubble time At = (t+ — t_) < H~^, 
gravity is negligible during this free fall. If we estimate the contribution of the gravity term 
in the Euler equation (|3.17|) during the transition we get (t+ — t^)dtip- = 0(At/cphys)a- = 
C»[(ifAt)(if/fcphys)]5- = 0[{HAt){H/kphys)]ip_ < ip- except for an initial phase v?- = 
mod TT (which leads to a dip in the spectrum). Thus, the fiuid is moving inertially in the 
sense of Newton, the velocity stays constant, and the density contrast grows linearly in time 

S = S_ + k{r]- , (3.30) 

where y/S-tp- = An sin (</?_) is the pecuhar velocity at 

Since we have no jumps in pressure po, the density contrast 6 and the fiuid velocity ip 
stay continuous throughout the whole transition, in particular at the matching points of the 
different regimes. Gravity remains negligible during the entire transition since the growth of 
Tp due to gravity at the end of the transition can be estimated by ip+—ip_ = C(Atfcphysa+) = 
0[HAt{H/kpi,ys)5+] = 0[{HAtytp_] <^ tp., in the limit under consideration. At the end 
of the transition this solution has to be joined to the pure radiation dominated regime for 
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T < Ti,. Since the amplitude of the density contrast grows hnearly during the transition, 
the final amplitude Aout = ^+ is enhanced linearly in modulated by the incoming phase 



A, 



out 



A 



m 




(3.31) 



with ki = \/3/Ar], At] = 77+ — ?]_. The envelope of the linearly growing peak structure 
for subhorizon scales starts at the scale ki which corresponds to a CDM mass of Mi = 
9 X IQ-^Mq. 



2. Lattice QCD 



For our fit to lattice QCD, Eq. (|2.^ ), the time evolution separates into four regimes. 
First we have a WKB regime, where the sound speed is slowly changing compared to the 
frequency of the oscillation, i.e. the WKB condition \dcs/ drj\/ Cg -C = kcg is satisfied. The 
WKB solution is 

cos ruj{r]')d7]' (3.32) 



(3c2)V4 J, 

VS^J = An(3c^)^/Sin r U{7]')d7]' 

Jo 



where Am, as before, denotes the incoming subhorizon amplitude in a radiation dominated 
universe. The decreasing sound speed in the WKB regime leads to an increase in the density 
contrast and a decrease in peculiar velocity compared to the purely radiation dominated case. 
Close to rj^ the WKB condition can no longer be satisfied: Since Cs decreases by some power 
in [i]- — 1]), Eq. ( p.25|) , the left hand side of the WKB condition diverges like l/(?7- — f]), 



whereas the right hand side for any fixed k goes to zero. Therefore, the WKB solution for 
a given mode is valid up to a certain time fj{k) < ?7_. r/(A;) approaches ?7_ for increasing k. 

Close to rj^ there is a second regime (overlapping with the WKB regime for large enough 
k), where Cg can be approximated by the power law of Eq. ( |2.25| ), 



rj_ — 7] 



Cs = C J (3.33) 

7 + 1 

where (5 is defined for convenience below. Inserting this power law into the oscillator equation 
( p.28| ) gives a Bessel differential equation. The new variable z is the argument of the cosine 



in the WKB expansion. 



V- 



z = k j' Cs{v')dv' = '2Cf3{kr]^) [^^^^^ j ■ (3-34) 

The index (3 of the Bessel functions is obvious from the ratio of the two fundamental solutions 
Jj^{z)/Nf^{z) oc z"^^ oc (?7_ — riY^f^ . Since the two fundamental solutions for — are 5 — >• 
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const and S oc (?7_ —t]), it follows that (3 = (3. The prefactor must be to make the second 
fundamental solution constant, 

5 = [-B^J(,{z) + B^Nf^iz)] . (3.35) 

The peculiar velocity follows via -ij) = S'/k, see Eq. ( ^■27| ). Matching the WKB solution 
with the Bessel solution in the overlap region z ^ 1, where Jpiz), Np{z) — > ^2/ (tiz) cos[z — 
(7r/2)(/3 + 1/2, 3/2)], gives the normalization B"^ = Bf + 5|, 

B = A,^^{2C(3)-^kv^y/'~f' . (3.36) 

At the beginning of the first-order transition, rj = •)]_, the Bessel solution (second regime) 
is matched to the third regime, i.e. the free fall regime with = 0. As in the bag model, 
peaks in the transfer function for large k are obtained for the mixing angles ip in the Bessel 
solution {Bi = Bcos(p,B2 = B simp) which give the maximal {ipl at the matching point. 
With Jp{z) {z/2)'^/T{l3 + 1) for ^ < 1 and Npiz) = [Jfsiz) cos(/57r) - J^p{z)]/ smilSn) we 
find 

^-'^ = An ^^..f.y . {kv-f-'^ . (3.37) 

The evolution in the third regime, = 0, is as before. The restoring force in the oscillator 
equation, Eq. ( p. 281) , vanishes, the fluid moves with constant peculiar velocity, ip{ri) = 
and the density contrast grows linearly in time. 

After the transition [rj > 77+, fourth regime) the universe is radiation dominated (c^ = 
1/3) and the solutions are oscillations with amplitude Aont = The amplification of 

peaks in the outgoing amplitude A^nt is given by ^outlpeaks ~ ip^^^k^rj, 



A\r 



{kr^^f-HA^. (3.38) 



peak. 31/4 r(/?+l)sin/5 



TT 



The dependence on the latent heat is encoded in the duration Arj of the transition. C and j3 
are determined by the normalization and the power law for 0^(77) shortly before the onset of 
the first-order phase transition, Eq. ( p.33| ). Replacing /5 by 7 (= the exponent in the QCD 
fit to the entropy density) the amplification factor is 



37+1 

A ■ ■ 



Aw 



peaks 




2(7+1) 

(3.39) 



for large k. Since < 7 < 1 the exponent of the power law varies between 1/2 and 1. The 
normalization of the peak amplification for large k is encoded in k2- In our fit to lattice 
QCD we take 7 = 1/3, hence f3 = 1/4 



A 



out 



A 



in 



peaks 




(3.40) 



For Rl = 0.2, k2 corresponds to a CDM mass of M2 = 2 X 1O-^°M0. The spectrum for this 
case together with the asymptotic envelope is plotted in Fig. 
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3. Crossover 



In the case of a crossover, the amplification is occurring for scales around the Hubble 
radius at the transition only. Subhorizon scales always stay in the WKB-regime and therefore 
the spectrum is fiat for these scales. However, the amplitude for subhorizon scales is damped 
during the phase transition. The same damping occurs in the case of a first-order phase 
transition. It has been neglected in the analytic discussion, since it is a small correction. 
The time evolution of the density contrast of subhorizon scales, including damping terms, 
reads 

8" + (1 + 3c^ - Qw)-5' + uj\r])8 = . (3.41) 
a 

The friction term in this equation can be eliminated if we write the density contrast as 
5 =: ^^^'^ u{ini). Inserting this ansatz into Eq. ( |3.41| ), using the Friedmann equation and 
w' = 3aH{l + w){w — Cg), we obtain an undamped oscillator equation with frequency 
07(77) = kcs{ri) for the quantity u{ri). For a crossover the sound velocity never vanishes, and 
subhorizon modes always satisfy the WKB-condition. Therefore, we obtain for the density 
contrast and the peculiar velocity 



^ = ~ (3c2)i/4 ^ ^ / ^(^)^^ ^3.42) 

a', 



o 1 /. Li- \/3(l + w) r 
^ = AUSc'.y/'^^^ ^sin / uiv)dv. 

The reduction of the amplitude of subhorizon scales is given by 

in 



An {Ha 




(3.43) 



Therefore the subhorizon amplitudes are reduced to 83 % of their initial value. This damping 
is a general feature appearing whenever the number of degrees of freedom is reduced. Also 
in the case of a first-order transition, this damping factor should be taken into account. 
However, it changes the huge amplifications merely by the constant factor 0.83. 



D. Analytic solution for CDM 

Cold dark matter (CDM) is assumed to be the major matter content of the universe 
today. At the time of the QCD transition, however, the contribution of CDM to the total 
energy density was negligible, pCDM^^RAD _ a/aequs,\ity ~ 10~^. We consider a type of CDM 
which is non-relativistic (p ^ p) at the QCD transition and which is only coupled via 
gravity to the radiation fiuid. CDM moves in the external gravitational potential provided 
by the dominant radiation fiuid. During the QCD transition, the big amplifications of the 
density contrast in the radiation fiuid S^^^ leads to a big amplification in the gravitational 
potential, i.e., the force term in the Euler equation for CDM. The CDM is accelerated to 
higher velocities at the end of the transition. Therefore, we also get peaks and dips in the 
cold dark matter fiuid. 
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The subhorizon evolution of CDM in a purely radiation dominated universe is just inertial 
motion. This can be seen from the Euler equation ( 3.17|) : The sound velocity of cold dark 



matter is zero, = 0, and the gravitational potential a is given using the Poisson equation 
( p.l8| ) for the dominant density contrast S^^^ Eq. ( p.30|) . In the leading order in x = k/Ti, 



the gravitational force can be neglected and the subhorizon evolution of the CDM velocity 
is obtained: 

^CDM ^ ^3 44^) 

X 

The velocity of CDM in a radiation dominated universe just redshifts to zero on subhorizon 
scales, see Fig. ^. C is an integration constant of order A[^. The corresponding evolution of 
the density contrast S'"^^ follows from the continuity equation (|3.16|) , 



rCDM 



C\nx + D. (3.45) 



This logarithmic growth of 5*^^^ can be seen in Fig. IC before and after the transition. The 
shape of S'"^^ can also be seen in the transfer functions Figs. ^^,0 on scales above the 
horizon scale Mi. 

During a first-order QCD transition the evolution of the gravitational potential is dom- 
inated by the linear growth in time of 5^"^^, Eq. ( |3.3CI| ). From the Poisson equation ( p.l8| ) 
follows that the gravitational potential grows linearly in time 

« = -^^(^-^-)(^)^?^°, (3.46) 

see Fig. p. As before, we will take the transition time to be short compared to the Hubble 
time At < H^^, hence Arj < IHr^. The gravitational force in the Euler equation ( |3 . 1 7| ) and 



hence the fluid acceleration for CDM also grows linearly in time during the transition. The 
fluid velocity grows quadratically 

^CDM ^ 3 [^^^]2 ^RAD ^ ^CDM_ ^3 47) 



The initial CDM velocity ip^^ is given by Eq. ( |3.44| ). On scales well below the horizon the 



velocity is redshifted to small values and can be neglected. However, on subhorizon scales 
closer to the horizon the initial velocity dominates, because [7i(A?7)]^ ^ 0.1(0.01) for the bag 
model (lattice fit) is small. According to the continuity equation (|3.16|) the density contrast 
^CDM g];Q-^ys cubically 

^CDM ^ ^CDM ^ 1 1^^^ [^^^]3 ^RAD_ ^3_4g) 

Due to the cube the growth of (JCDm (j^ri^g the QCD transition is small. 

The major amplification effect comes from the higher velocity at the end of the transition 
(see Fig. |^), which leads to 

<5CDM^^^ ^ ^™/tr7+ In(^) + 5™. (3.49) 
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The amplification of tlie density contrast during tlie transition is negligible compared to the 
enhancement of the velocity. This velocity enhancement during the transition leads to an 
additional logarithmic growth of 5^^^ after the transition, 



rCDM 



iv) 



^ f^v\ .7, RAD 



4 V r/^ 



A:r^+ln(^) + 5™. 



(3.50) 



Here we have used ( p.47| ) to obtain the leading subhorizon approximation. This additional 
logarithmic growth of 5^^^ is most clearly seen in Fig. 

The number of CDM peaks in the numerical result, Figs. |^|^, is only half of the radiation 
peaks. This comes from the fact that the peaks in (JCDm jjg logarithmic curve, 

which is still the main contribution at the scales plotted. The dips correspond to modes, 
with a maximal negative amplification. 

The jump in the gravitational potential at the end of the QCD transition in Fig. ^follows 
from the Poisson equation ( p.l8|) . Since 6 is continuous, but jumps, the right hand side 
abruptly changes by a factor 2. Therefore a jumps by the same factor. 

Let us compare the 'QCD peaks' in the CDM spectrum with the CDM spectrum without 
phase transition. In the limit pCDM^^RAD ^ ^^iq exact solution to Eqs. (|3.1CI|) - (|3.12|) in 
the radiation dominated universe without transition reads {x = k/Ti): 



rCDM _ 3^in 



CDM 



— sinC^) + cos(^) 



x2 + 6lv^ 'V3 



1 - 



6 



, sin(— ^) + cos(— ^) 



In the subhorizon limit this solution reduces to 



(3.51) 
(3.52) 



rCDM 



CDM 



3v4ir 
3Ai„ 1 



2 X 



30.077 



which fixes the constants C and D in Eqs. ( p.44|) and (|3.45| ). 
We are now able to calculate the enhancement factor 



5. 



transition 



5, 



no transition 



CDM 



(3.53) 
(3.54) 



(3.55) 



Let us evaluate S at equality, i.e., rj^Jri^ ^ 2 x 10^: 

^(r?cq) ~ 1 + 



1 /A?;' 



2 



RAD 



fcr,X-)/ln(-^, . 



(3.56) 
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The enhancement at the smallest scales M'~^^^ ~ 10 ^°M0 to which our calculation applies 
and at some intermediate scale M*^^^ ~ 10~^^Mq reads 

^(-J - { l:?,rm ' } io--(io-'')M. . (3.57) 

where we have used Eq. ( p.37| ) for the case of the lattice fit. This shows that both models 
lead to large enhancements of the CDM density contrast at equality for small enough scales. 
For the lattice fit the enhancement becomes important at scales below M^^^ ~ 10~^^Mq. 



IV. IMPLICATIONS CLUMPS IN CDM 



A. Fate of the peaks in the radiation fluid 



Before discussing CDM, let us show that the large peaks in the density spectrum of the 
hadron-lepton-photon fiuid do not lead to any observable consequences. 



1. CoUisional Damping at Neutrino Decoupling 



The acoustic oscillations in the radiation fiuid get damped by neutrino diffusion at the 
time of neutrino decoupling. This damping is analogous to Silk damping at photon de- 
coupling. The muon resp. tau neutrinos, which are coupled to the RAD fluid via neutral 
current interactions only, decouple at Tj^'^'^,^ ~ 2.2 MeV from the Hubble scale Rn, which 
follows from Ref. [0. The electron neutrinos interact by charged and neutral currents and 
decouple slightly later, T^^'^ ~ 1.4 MeV. By the time of neutrino decoupling at the Hubble 
scale all inhomogeneities in the RAD fluid on scales below ~ 10~^Mq in CDM mass are 
wiped out by neutrino diffusion (cf. the QCD horizon scale is 10~^Mq in CDM mass), as 
shown below, Eq. ( [4.5| ). Therefore our QCD peaks cannot affect big bang nucleosynthesis 
(BBN). 

It is important to distinguish the total decoupling of neutrinos, i.e. neutrino decoupling 
at the Hubble scale, when the neutrino scattering rate F is less than the Hubble rate, 
Tu/H < 1, from the neutrino decoupling with respect to a certain mode given by Fj^/wphys < 
1, when the typical neutrino scatters less than once during an acoustic oscillation time of 
one particular mode. The mode- dependent decoupling temperature T^^^{k) is related to the 
total decoupling temperature by 

This follows because the neutrino interaction rates Fj, are proportional to and fcphys c< 
1/a (xT, hence Ft./cijphys oc T^. 

To compute the damping of acoustic oscillations in the radiation fluid by neutrino diffu- 



sion we follow Weinberg |47]. For a radiation fluid shear viscosity is dominant, bulk viscosity 



vanishes and heat conduction is negligible. The shear viscosity is given by 
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4 

^visc = Y^$]p^,i-^, • (4.2) 

p^^ denotes the energy density of a neutrino species, t^. is the typical colhsion time. In the 
subhorizon hmit the Navier-Stokes equation and the continuity equation give 

S" + ^r^^.^^k5' + u;H = Q, (4.3) 

Ptot 

a damped oscillator. The damping factor for the mode at a given conformal time t] is 

X /■''?max 

D{k, if) = exp [-- {kpi,ys/ptot)Vvisckdr]'], (4.4) 

Here the upper limit of the integral is the conformal time rjmax = Min [ri,ridec{k)], because 
collisional damping of the 7-^^-hadron fluid by neutrinos ceases at decoupling of the mode 
k considered. The damping per oscillation is largest for ujt,^ = uo/Ty ^ 1, i.e. immedi- 
ately before neutrino decoupling for a given mode. But note that subhorizon modes get 
strongly damped long before the mode decouples from neutrinos, because a weak damping 
per oscillation is compensated by many oscillations per Hubble time. 

For a first application we ask what modes are already damped by the time of the QCD 
transition. At the QCD transition, T = T^, the interaction rates for electron and muon 
neutrinos (resp. antineutrinos) with the leptons are the same, T^^ = Tp^ = T^^ = Tp^ = 
S.lGpT^ |^6|, since electrons and muons are still relativistic. The r— neutrinos interact only 
via neutral currents with the leptons and have a lower interaction rate, T,^^ = Tp^ = O.GGpT^. 



Evaluating the damping integral, Eq. (4^), at T = we find that the damping factor 
D{k,Ti,) is < 1/e for (/cphys/-ff)r* > 10"^, which corresponds to Mcdm < lO'^^M©, i.e. 
acoustic oscillations on these small scales are wiped out before the QCD transition. Therefore 
no peaks in the RAD or in the CDM transfer function can develop below this scale. This 
small-scale cut-off is independent of the bubble separation scale. 

Next we consider T < T^^^[H) f» 1.4 MeV and evaluate the final damping factor D{k). 
At the time of neutrino decoupling we take a purely radiation dominated universe, consisting 
of 7, and v' s only to evaluate the damping. Muons have disappeared since ^ T and 
since ^ tn- The interaction rate for electrons is still given by charged and neutral 
currents, T^^ = Tp^ = l.SGpT^, whereas muon and tau neutrinos have the same lower 
(neutral current) rate, T^^ = Tp^ = T^^ = Tp^ = 0.3GpT^ [^. The final damping of a 



certain scale (expressed in invariant CDM mass) from Eg. ( [4.4|) is 

^ = exp[-(^) ], (4.5) 

with the neutrino damping scale Md ^ 1.9 x W'^Mq in CDM mass. This is 3 x 10~^M^^^ 
at T = 1.4 MeV and corresponds to length scales A = 1/30-Rh- 

We conclude that the large peaks in the spectrum of the radiation fluid are damped away 
by the time of BBN. The large peaks in the RAD fluid generated by the QCD transition do 
not give a new mechanism for inhomogeneous BBN. 
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2. Black Hole Formation at the QCD Transition? 



It was suggested in the literature [ PU| , P^^ that the QCD transition could lead to forma- 
tion of IMq black holes, which could account for dark matter today. Jedamzik ||T9| proposed 
to identify such black holes with the MACHOs observed by microlensing [^]. He pointed 
out that the formation of black holes should be particularly efficient during the QCD epoch 
due to the significant decrease in the effective sound speed. 

In order to form a black hole in a radiation dominated universe, the density contrast 
inside the Hubble radius should be in the range 1/3 < 5^ < 1 W^- For an observable 



< 5h < 1 _ 

amount of IM© black holes today, i.e. fign = ^(1); the fraction of energy density converted 
to black holes at the QCD transition must be C(aQCD/aequaiity) ~ 10~®. For a Gaussian 
distribution this requires (5rms ~ 0.06 (without including any enhancement from the QCD 
transition) |5^. The QCD transition gives an enhancement factor (at the horizon scale) 
of 2 for the bag model and of 1.5 for lattice QCD in our linear perturbation treatment. 
Figs. § - 0. This indicates a corresponding reduction in the required preexisting perturbation 
spectrum at the solar mass scale. Cardall and Fuller used a qualitative argument of Carr 
and Hawking |53| and the bag model and also obtained a factor 2 reduction in the required 
preexisting perturbation spectrum [Q]. These QCD factors of 1.5 or 2 are so modest that a 
preexisting Harrison-Zel'dovich spectrum with COBE normalization is very far from giving 
a cosmologically relevant amount of black holes |2^. One would have to put in a fine-tuned 
tilt {n — 1) 0.4 to get the desired amount of black holes. However, this tilted spectrum 
would overproduce primordial black holes on scales which are only a factor 50 below the 
Hubble radius at the QCD transition. Therefore a break in the preexisting spectrum below 
the QCD scale would be required, a second fine tuning. 

We conclude that the QCD transition enhances the probability of black hole formation, 
but the preexisting spectrum needs to be fine tuned around the QCD scale and the major 
effect would not be due to the QCD transition. 



B. Kinetic decoupling of CDM 



In Sec. HI we established the generation of peaks in the CDM spectrum during a first- 
order QCD transition. Below, we discuss the properties of the most prominent CDM candi- 
dates at temperatures of the QCD scale. We show that weakly interacting massive particles, 
like the lightest supersymmetric particle or a heavy neutrino, are kinetically tightly coupled 
to the radiation fiuid at the QCD transition and are included in the radiation fiuid in our 
equations. Two CDM candidates that are kinetically decoupled from the radiation fiuid and 
have p = at the QCD epoch are the axion and primordial black holes. In our figures and 
equations kinetically decoupled CDM have been labeled by CDM for brevity. 



1. Neutralinos and heavy neutrinos 

We start our discussion with the lightest supersymmetric particle In the minimal 

supersymmetric standard model this is the neutralino (we assume it is stable). Constraints 
from LEP 2 and cosmology, together with the assumption of universality at the GUT scale. 
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show that its mass is > 42 GeV It is essential to distinguish between the chemical 
freeze-out and the kinetic decoupling of neutralinos. The chemical freeze-out determines the 
amount of neutralinos today, and it happens when the annihilation rate of neutralinos drops 
below the Hubble rate, FannZ-f^ < 1- When the neutralinos become nonrelativistic the rate 
for neutralino annihilation, Fann = {vcrann)n^, is suppressed by the Boltzmann factor in the 
number density of the neutralinos, ~ (m^T)^/^ exp(— m^^/T). The freeze out temperature 
of the neutralino |^ is Tf ~ > 2 GeV, and neutralinos are chemically decoupled at 

the QCD transition. 

Kinetic decoupling, in contrast, is determined by the elastic scattering between neu- 
tralinos and the dominant radiation fluid. The interaction rate for elastic scattering is 
Fei = r^J^i = {vaci)n, where n ~ is the number density of relativistic particles, e.g. elec- 
trons or neutrinos. An order of magnitude estimate shows that aei is similar to the cross 
section for the elastic scattering of neutrinos and neutrons, Uei ~ GpT^ We have to 



distinguish between the regime of perfect kinetic coupling, i.e. neutralinos tightly coupled 
to the radiation fluid, an intermediate regime where the neutralinos scatter elastically but 
the number of collisions is not sufficient to drag them along dissipationless, and the regime 
of kinetic decoupling, T^i/H < 1, which is roughly at T < 1 MeV, since the neutralino 
interacts weakly. 

Let us estimate the regime where the neutralinos belong to the perfect radiation fluid. 
Perfectness of a fluid (dissipationless) always refers to an external time scale, uj~^ of an 
acoustic oscillation or the Hubble time, etc. A fluid behaves as a perfect fluid, if t^Trdax 1, 
i.e. the external time scale is larger than the relevant relaxation time, the fluid is contin- 
ually in local thermal equilibrium, and no entropy production occurs. For the coupling of 
neutralinos to the radiation fluid the relaxation time is given by Treiax = A^TcoU, where is 
the number of collisions needed to completely change the momentum of the neutralino due 
to collisions with the radiation fluid. The momentum transfer at a collision of a lepton with 
the neutralino is of order p\ ~ T. The kinetic energy of the neutralino is given by equipar- 



tition, p^/2m^ ~ T, hence its momentum is ~ ^^m^T. The fractional change of the 

neutralino momentum from one collision at the QCD transition is ^"PxlVx ~ \Jt /fn^ <^ 1. 

After A^ collisions the total rms change of momentum is ('^Px/^'x)rms ~ y^A^T/m^. Local 
thermal equilibrium is obtained if the bulk motion of the neutralinos is governed by the 
leptons, i.e. the fractional change of the neutralino momentum is of order 1. The number 
of collisions needed to completely change direction is A^ ~ '^x/-^ ~ "^^^ ~ GeV 

and T = T^. The collision time is given by the weak interactions rate F„ = r^^^y ~ GpT^, 
and the relaxation time is given by 

W = ATreou ~ AT X (4.6) 

with A^ ~ m^/T. If we compare the relaxation time with the frequency of the acoustic 
oscillations, we find that the condition for a perfect fiuid, tuTreiax <^ 1, at the QCD transition 
is satisfied for scales A > Adec(^*) ~ 10~'^i?H ('^x ~ ^0 GeV). Hence the neutralinos on these 
scales are part of the radiation fiuid at the QCD transition. Below this scale, the neutralinos 
cannot follow the acoustic oscillations. On the other hand on the Hubble scale the perfect 
kinetic coupling of neutralinos to the radiation fiuid stops when the required relaxation time 
becomes more than a Hubble time, Trdax > ^h- This gives a temperature of T^'^"^ ~ 10 MeV. 
Down to this temperature neutralinos on the Hubble scale belong to the radiation fiuid. 
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Another CDM candidate is a heavy neutrino with mz /2 < mu < 1 Te V , where only 
the upper mass bound gives cosmologically relevant CDM. Again, the kinetic decoupling of 
these neutrinos happens way below T^. Therefore, we find that the neutralino or a heavy 
neutrino would be tightly coupled to the radiation fluid during the QCD transition. 

2. Axions 



One candidate for CDM at the QCD epoch is the axion [jS^]. We therefore find a new 
mechanism to make axion clumps. 

Axions could be the dominant matter today if their mass is small, i.e. ma ~ 10"^ eV, 
which corresponds to a breaking of the Peccei-Quinn (PQ) symmetry at the scale /pq ~ 10^^ 



GeV These axions could be produced coherently due to the initial misalignment of the 
axion field and by the decay of axionic strings. The initially misaligned axion field starts 
to oscillate coherently when the axion mass has grown to ma{Ti) ~ 3H{Ti), where Ti ~ 1 



GeV |58|. Thus, below Ti ~ 1 GeV the oscillating axion field evolves as CDM, pa oc a~ . 

In a first-order QCD transition axionic CDM falls freely into the large potential wells 
from the dominant radiation fluid. This produces peaks and dips in the spectrum of axionic 
CDM. These peaks trigger the formation of clumps of axions with masses 10~^° - 10~^°Mq 
shortly after equality, see Sec. IV. C. 

Another mechanism to produce clumps of axions has been suggested in Ref. [^. If 



the reheating temperature after inflation is above /pq, the axion field is inhomogeneous on 
scales larger than Ri{{Ti). Inhomogeneities in axion density evolve into axion miniclusters 
of typical mass Mmc ~ IO^^^Mq and typical radius R^c ~ 10^ m ~ O.IRq ~ 10"^ AU. 
These axion miniclusters may be enhanced by the free-fall of matter during a first-order 
QCD transition. For these non-linear axion inhomogeneities a quantitative result cannot be 



obtained from the linear analysis of Sec. III. Let us emphasize that, in contrast to Ref. ^9 
we predict axion clumps independently from the reheating temperature after inflation. 

3. Primordial black holes formed before the QCD transition 

A further possibility for CDM that decoupled kinetically long before the QCD transition 
is primordial black holes (PBH) produced before the QCD transition and therefore with 
masses Mbh IMq. In order to survive until today PBH should have Mbh > lO^^g 
pa 10~^^Mq. PBH in the range from IO'^^Mq to IO'^^Mq would radiate too strongly to 



be compatible with 7-ray observations ||6^ . The production of PBHs arises naturally from 
hybrid inflationary scenarios [|6l| ,^. 

Our analysis of the QCD phase transition predicts that these small black holes would 
be attracted by the large potential wells in the radiation fluid and would therefore tend to 
form miniclusters of masses up to 10~^^ Mq. 



C. Clumps in CDM 

CDM in the form of axions or PBHs is not subject to damping as the radiation fluid, 
thus the peaks in CDM will survive until structure formation starts. The free streaming 
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scale of CDM is way below our smallest scales, because the initial velocity of axions or PBHs 
is completely negligible. An overdensity of CDM in the form of axions or PBH (or any other 
matter that is kinetically decoupled at the QCD epoch) decouples from the cosmic expansion 
when its density contrast becomes nonlinear, {5p/ p)r ~ 1 (condition for turn around). It 
collapses and virializes by violent gravitational relaxations and forms a clump of CDM. For 
spherical collapse the final viral radius is half of the radius at turn around [ |6^ . 

We take a COBE |^ normalized spectrum and allow for a tilt |^ — 1| < 0.2. During the 
radiation dominated regime, {5p/ p)k for CDM continues to grow logarithmically. In {5p/ p)\ 
another logarithm comes from summing modes up to A; ~ l/i?, where R is the radius of 
the window function. The enhancement factor £ of CDM density fluctuations during the 
QCD transition has been obtained in Eq. (|3.55| ). Putting these factors together we obtain 
the amplitude of CDM perturbations of size R at the time of equality: 



CDM 



n-1 
2 



1 ( ^ 



^fc(teq) (4.7) 





/co denotes the comoving wave number of the mode crossing the horizon today. In the 
following we assume Vt^h'^ = 1/4, and thus ~ 6000. For a CDM mass of 1O~^^(1O~^°)M0 
a tilt of 72 — 1 = 0.2 gives a factor ^ 16(23). The logarithms contribute another factor 
^ 94(140). The enhancement factor has been calculated in Sec. Ill to be S{teq) ~ 2(8) for 
the lattice QCD fit. Looking at 3a peaks we find without tilt (n — 1 = 0) 

3(T,n=l 

(teq) ^ 0.1(0.6) for Md,^p = 1O-1^(1O-2°)M0 . (4.8) 

R 

This implies that these clumps go nonlinear at z^i ~ 600(3600) and collapse to clumps of 
radius i?ciump ~ 14(O.O4)i?0. With tilt (n - 1 = 0.2) we find 

3cr,n=1.2 

(teq) ^ 2(16) for Mci,^p = 1O-1^(1O-2°)M0 . (4.9) 

R 

These over-densities start to collapse even before teq because locally CDM starts to dominate 
at ~ 2(16)Teq. This leads to clump radii of -Rdump ~ 0.7(0.003)/?©. 

We conclude that the peaks in the CDM spectrum lead to clumps of masses 10~^° — 
10^^'^ M0. Today, these clumps would have a density contrast of 10^° — 10^^, where the lower 
value corresponds to a 10~^^Mq clump from an untilted CDM spectrum, the bigger value 
is for a 10~'^^Mq clump from a tilted CDM spectrum. The evolution of these clumps in the 
late stages of structure formation remains to be investigated (disruption, mergers, etc.). 

For larger enhancement, e.g., if it should turn out that the latent heat is bigger than the 
value from present lattice QCD calculations, more compact clumps are possible. These could 
be subject to femto-lensing With the values of our lattice fit, the CDM clumps are 

not compact enough to lie within the Einstein radius, which is -Re ~ 0.02i?Q for a 1O~^^M0 
clump. 

The clumping of CDM changes the expected reaction rates for some dark matter searches, 
because some of the rates depend on the space-time position of the detector, star, or planet. 
Especially experiments looking for axion decay in strong magnetic fields [^,0] would not 
yield limits on the axion mass. Maybe these experiments just tell us that we are not sitting 
in an axion cloud currently. These consequences remain to be studied further. 



30 



ACKNOWLEDGMENTS 



We are very grateful to U. Heller and F. Karsch for providing us with their lattice QCD 
data sets. We like to thank J. A. Bardeen, V. Berezinsky, B. J. Carr, U. Heller, J. Ignatius, 
P. Jetzer, F. Karsch, H. Markum, V. Mukhanov, J. Silk, P. Sikivie, and N. Straumann for 
helpful discussions and references to the literature. C.S. thanks J. Silk and the Center for 
Particle Astrophysics in Berkeley for hospitality. D.J.S. and P.W. thank the Swiss National 
Science Foundation for financial support. D.J.S. is an Alexander von Humboldt fellow. 



31 



REFERENCES 



[1] For recent reviews see K. Kanaya, Nucl. Phys. B (Proc. Suppl.) 47, 144 (1996); E. 

Laermann, Nucl. Phys. B (Proc. Suppl.) 63, 141 (1998). 
[2] Y. Iwasaki et al., Z. Phys. C 71, 343 (1996); Nucl. Phys. B (Proc. Suppl.) 47, 515 

(1996). 

[3] Y. Iwasaki et al., Phys. Rev. D 46, 4657 (1992); 49, 3540 (1994); B. Beinhch, F. Karsch, 

and A. Peikert, Phys. Lett. B 390, 268 (1997). 
[4] C. Bernard et al., Phys. Rev. D 54, 4585 (1996). 
[5] J. W. Harris and B. Miiller, Ann. Rev. Nucl. Part. Sci. 46, 71 (1996). 
[6] J. Ignatius et al, Phys. Rev. D 49, 3854 (1994); 50, 3738 (1994). 
[7] M. B. Christiansen and J. Madsen, Phys. Rev. D 53, 5446 (1996). 
[8] E. Witten, Phys. Rev. D 30, 272 (1984). 

[9] K. Sumiyoshi and T. Kajino, Nucl. Phys. B (Proc. Suppl.) 24, 80 (1991); P. Bhattachar- 

jee et al., Phys. Rev. D 48, 4630 (1993). 
[10] J. H. Applegate and C. J. Hogan, Phys. Rev. D 31, 3037 (1985); J. H. Applegate, C. J. 

Hogan, and R. J. Scherrer, Phys. Rev. D 35, 1151 (1987); G. M. Fuller, G. J. Mathews, 

and C. R. Alcock, Phys. Rev. D 37, 1380 (1988); R. A. Malaney and G. J. Mathews, 

Phys. Rep. 229, 145 (1993). 
[11] In-Saeng Suh and G. J. Mathews, e-print |astro-ph/9805T79| (1998). 
[12] C. Schmid, D. J. Schwarz, and P. Widerin, Phys. Rev. Lett. 78, 5468 (1997). 
[13] E. R. Harrison, Phys. Rev. D 1, 2726 (1970); Ya. B. Zel'dovich, Mon. Not. R. Astr. Soc. 

160 IP (1972). 

[14] V. Mukhanov and G. Chibisov, Pis'ma Zh. Eksp. Teor. Fiz. 33, 549 (1981) [JETP Lett. 
33, 532 (1981)]; A. Starobinsky, Phys. Lett. B 117, 175 (1982); A. Guth and S.-Y. Pi, 
Phys. Rev. Lett. 49, 1110 (1982); S. Hawking, Phys. Lett. B 115, 295 (1982). 

[15] J. M. Bardeen, P. J. Steinhardt, and M. S. Turner, Phys. Rev. D 28, 679 (1983). 

[16] J. M. Bardeen, in Cosmology and Particle Physics, ed. A. Zee (Gordon and Breach, 
New York, 1989); see also J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[17] V. F. Mukhanov, H. A. Feldman, and R. H. Brandenberger, Phys. Rep. 215, 203 (1992). 

[18] G. F. Smoot et al., Astrophys. J. 396, LI (1992). 

[19] K. Jedamzik, Phys. Rev. D 55, 5871 (1997). 

[20] G. Jungman, M. Kamionkowski, and K. Griest, Phys. Rep. 267, 195 (1996). 
[21] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983); C. Hagmann et al., Phys. Rev. Lett. 80, 
2043 (1998). 

[22] MACHO Collaboration, Astrophys. J. 486, 697 (1997); EROS Collaboration, Astron. 

Astrophys. 324, L69 (1997). 
[23] P. Widerin and C. Schmid, e-print |astro-ph/980TT63| (1998). 
[24] D. J. Schwarz, preprint |gr-qc/9709027| (1997). 

[25] A. Chodos et al., Phys. Rev. D 9, 3471 (1974); T. De Grand et al., ibid. 12, 2060 (1975); 
T. D. Lee, Particle Physics and Introduction to Field Theory (Harwood Academic Pub- 
lishers, Chur, 1981). 

[26] G. Boyd et al., Phys. Rev. Lett. 75, 4169 (1995); Nucl. Phys. B469, 419 (1996). 

[27] MILC Collaboration, Phys. Rev. D 55, 6861 (1997). 

[28] J. P. Blaizot and J. Y. Ollitrauh, Phys. Rev. D 36, 916 (1987). 

[29] J. Engels et al., Phys. Lett. B 396, 210 (1996). 



32 



[30] E. Farhi and R. L. Jaffe, Phys. Rev. D 30, 2379 (1984). 

[31] R. Pisarski and F. Wilczek, Phys. Rev. D 29, 338 (1984); F. Wilczek, Int. J. Mod. Phys. 

A 7, 3911 (1992); K. Rajagopal and F. Wilczek, NucL Phys. B399, 395 (1993). 
[32] F. R. Brown et ah, Phys. Rev. Lett. 20, 2491 (1990). 
[33] See e.g. D. H. Rischke and M. Gyulassy, Nucl. Phys. A 597, 701 (1996). 
[34] T. DeGrand and K. Kajantie, Phys. Lett. 147B, 273 (1984). 
[35] G. M. Fuller, G. J. Mathews, and C. R. Alcock, Phys. Rev. D 37, 1380 (1988). 
[36] L. D. Landau, E. M. Lifschitz, and L. P. Pitajewski, Lehrbuch der theoretischen Physik, 

vol. V, 7. ed., Berlin 1987. 
[37] L. P. Csernai and J. I. Kapusta, Phys. Rev. D 46, 1379 (1992). 
[38] M. Hackel et al., Phys. Rev. D 46, 5648 (1992). 
[39] H. Kurki-Suonio, Nucl. Phys. B255, 231 (1985). 
[40] K. Kajantie and H. Kurki-Suonio, Phys. Rev. D 34, 1719 (1986). 
[41] K. Kajantie, Phys. Lett. B 285, 331 (1992). 

[42] M. Nagasawa and J. Yokoyama, Prog. Theor. Phys. 97 (Lett.), 173 (1997). 

[43] M. Nagasawa and J. Yokoyama, e-print |astro-ph/9612014| v2 (1997). 

[44] J. Hwang, Astroph. J. 415, 486 (1993). 

[45] C.L. Bennett et al, Astroph. J. 464, L1-L4 (1996). 

[46] A. Heckler and C. J. Hogan, Phys. Rev. D 47, 4256 (1993). 

[47] S. Weinberg, Astrophys. J. 168, 175 (1971). 

[48] D. L. Tubbs and D. N. Schramm, ApJ. 201, 467 (1975). 

[49] M. Crawford and D. N. Schramm, Nature 298, 538 (1982). 

[50] K. Jedamzik, e-print |astro-ph/980"5Ti7| (1998). 

[51] D. K. Nadezhin, I. D. Novikov, and A. G. Polnarev, Astron. Zh. 55, 216 (1978) [Sov. 
Astron. 22, 129 (1978)]; G. V. Bicknell and R. N. Henriksen, Astrophys. J. 232, 670 
(1979). 

[52] J. S. Bullock and J. R. Primack, Phys. Rev. D 55, 7423 (1997). 

[53] B. J. Carr and S. W. Hawking, Mon. Not. R. Ast. Soc. 1 68, 399 (1974). 

[54] C.Y. Cardall and G. M. Fuller, e-print |astro-ph/980TT0^ (1998). 

[55] J. Elhs et al., Phys. Lett. B 413, 355 (1997); e-print |hep-ph/9801ii5| (1998). 

[56] J. Ellis, in: Cosmology and Large Scale Structure, Les Houches, session LX, eds. R. 

Schaeffer et al. (Elsevier Sci., Amsterdam, 1996) pp. 825. 
[57] J. E. Kim, Phys. Rep. 150, 1 (1987); M. S. Turner, Phys. Rep. 197, 67 (1990); G. G. 

Raffelt, Phys. Rep. 198, 1 (1990); preprint |astro-ph/ 97072^ (1997); P. Sikivie, preprint 

hep-ph/ 96 113391 (1996). 
[58] D. Gross, R. Pisarski, L. Yaffe, Rev. Mod. Phys. 53, 43 (1981); M. S. Turner, Phys. 

Rev. D 33, 889 (1986). 
[59] C. J. Hogan and M. J. Rees, Phys. Lett. B205, 228 (1988); E. Kolb and I. I. Tkachev, 

Phys. Rev. Lett. 71, 3051 (1993); Phys. Rev. D 49, 5040 (1994); Astrophys. J. 460, 

L25 (1996). 

[60] B. J. Carr and J. E. Lidsey, Phys. Rev. D 48, 543 (1993); B. J. Carr, J. Gilbert, and J. 

Lidsey, Phys. Rev. D 50, 4853 (1994). 
[61] L. Randall, M. Soljacic, and A. H. Guth, Nucl. Phys. B 472, 377 (1996); preprint |Eep 



ph/9601296| (1996); J. Garcia-Belhdo, A. Linde, and D. Wands, Phys. Rev.D 54, 6040 
(1996). 



33 



[62] T. Padmanabhan, Structure Formation in the Universe (Cambridge University Press, 

Cambridge, 1993), chapter 8. 
[63] A. Gould, Astrophys. J. 386, L5 (1992); A. Ulmer and J. Goodman, ibid 442, 67 (1995). 



34 



FIGURES 



1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 





- 


1 


1 

P/Pideal 


1 

Nf 


1 

= 0,T^ 


1 

= 140 MeV 1 1 

= 260 MeV - 

bag model ■ ■ ■ ■ 






- 
















P/Pideal 










^ — 

1 1 1 1 1 



1.5 



2.5 



3.5 



FIG. 1. The energy density and the pressure of hot QCD relative to the energy density of an 
ideal quark-gluon plasma are plotted function of T/T^. The lattice QCD data for iV/ = ^ 
and Nf = 2 QCD Q are compared with the predictions of the bag model. 
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FIG. 2. The entropy density of hot QCD relative to the entropy density of an ideal QGP is 
shown for the same lattice QCD results as in Fig. |^. We also plot the entropy density for our fit 
to the lattice QCD data, for the bag model, and for a smooth crossover. 



35 



1 1 


1 1 1 1 1 1 

bag model ■ ■ ■ ■ 




lattice fit 


- 


crossover 










1 1 


1 1 1 1 1 1 



0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

a/a_|. 



FIG. 3. The behavior of w = p/p during the QCD transition as a function of the scale factor 
a. Although the pressure is reduced, it stays positive throughout the transition. 
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FIG. 4. The behavior of the sound speed = {dp/dp)s during the QCD transition as a 
function of a. For a first-order transition (lattice fit and bag model) the sound speed vanishes. 
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FIG. 5. The modifications of the density contrast for kinetically decoupled CDM (Hke ax- 
ions or primordial black holes) A'~"^^ = l^*""^^! (T^^/IO) and of the radiation fluid amplitude 
^RAD ^ -|- '^^'kad)^^'^ to QCD transition in the bag model. Both quantities are 

normalized to the pure Harrison-Zel'dovich radiation amplitude. On the horizontal axis the wave 
number k is represented by the CDM mass contained in a sphere of radius n/k. 
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FIG. 6. The same as Fig. ^, but for our lattice QCD fit. The straight line denotes the 
asymptotic envelope for small scales oc k^^^. 
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FIG. 7. The same as Fig. but for a QCD crossover. 
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FIG. 8. The time evolution of the density contrast, (5^"^^, and the pecuhar velocity, ji/j^^^ , of 
the radiation fluid for the mode k/ki = 7 in uniform expansion gauge. During the QCD transition 
in the bag model — marked by the 2 vertical lines — the velocity stays approximately constant and 
the density contrast grows linearly. The amplitude is normalized to 1 long before the transition. 
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FIG. 9. The time evolution of the peculiar velocity, ^/^ , of the kinetically decoupled 
cold dark matter fluid (axions or primordial black holes) and the external gravitational potential 

k 

gradient -^jp-oi for the mode k/ki = 37 in uniform expansion gauge, a is provided by the radiation 
fluid. During the QCD transition in the bag model — marked by the 2 vertical lines — the cold 
dark matter fluid is accelerated by the gravitational force. 
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FIG. 10. The time evolution of the CDM density contrast, 5^^^, corresponding to the CDM 
velocity of Fig. |9[ The major amplification of S^^"^ is due to the higher peculiar velocity at the 
end of the transition. 



39 



